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Abstract 

We give a nonrigorous derivation of the nonlinear Boltzmann equa- 
tion from the Schrodinger evolution of interacting fermions. The argu- 
ment is based mainly on the assumption that a quasifree initial state 
satisfies a property called restricted quasifreeness in the weak cou- 
pling limit at any later time. By definition, a state is called restricted 
quasifree if the four-point and the eight-point functions of the state 
factorize in the same manner as in a quasifree state. 

1 Introduction 

The fundamental equation governing many-body quantum dynamics, the 
Schrodinger equation, is a time reversible hyperbolic equation. The quantum 
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dynamics of many-body systems, however, are often modelled by a nonlinear, 
time irreversible (quantum) Boltzmann equation, which exhibits a particle- 
like behavior. This apparent contradiction has attracted a lot of attention 
over the years and, in particular, Hugenholtz [7] gave a derivation based 
on a perturbation expansion involving multiple commutators. He selected a 
class of terms from this expansion and argued that it gives the Boltzmann 
equation. To the second order in the coupling constant, Hugenholtz's claim 
was proved by Ho and Landau jHj. Beyond that it is not even clear that 
these terms satisfy the Boltzmann equation order by order, partly due to the 
complicated selection rules. 

In this paper, we present a derivation of the quantum Boltzmann equation 
under the main assumption that in the weak coupling limit (jSJ) the four- 
point and the eight-point functions of the state factorize at any time in the 
same manner as in a quasifree state, see (}2*2*j) and (|23|). A state with such 
a factorization property is called a restricted quasifree state. To rigorously 
verify this assumption, one has to analyze the connected m-point functions, 
a very difficult problem in our view. So it might appear that we have not 
improved much beyond the work [7] . Our approach however has the following 
two main merits. First: It identifies the concept, the restricted quasifreeness, 
to replace the independence in the classical setting so that the structure 
of the collision term, i.e., the quartic nonlinearity and the product of the 
factors F and 1 — F in the equation (JTUJ), appears as a simple consequence 
of this assumption. Second: Unlike Hugenholtz's approach which is tied to 
the commutator expansion, the restricted quasifreeness can now be verified 
using other methods such as field-theoretical techniques. 

Recent work by Benedetto, Castella, Esposito and Pulvirenti [2] has given 
an interesting different derivation. We had learned this work in a recent 
meeting and had subsequently sent them an early version of this manuscript. 

2 Definitions of the Dynamics 

We describe the quantum dynamics in the second-quantized formulation. For 
definiteness, we shall restrict ourselves to a fermion system. Our derivation 
is valid for bosons as well with the only difference being that some ± signs 
change along the derivation and in the final quantum Boltzmann equation 
(the terms (1 — F) change to (1 + F) in (JTUJ) ) . It should be noted that 
many-boson systems are in general more difficult to control rigorously than 
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many-fermion systems. The quantum Boltzmann equation for fermions also 
preserves the property < F < 1. On the contrary, the equation for bosons 
may blow up in finite time. 

Most of our setup is fairly standard and we recall the details briefly 
here. For background, see E] • The configuration space is a discrete torus 
A = Z d /LZ d , of a very large sidelength, LgN, which is kept finite through- 
out the argument. The Hilbert space for the fermions is the standard Fock 
space JF A := ® n>0 /\""Ha, where TC\ '■= £ 2 (A,C). 1 Because H\ is finite- 
dimensional, the same holds for T\. 

We shall work in momentum space, which for our finite lattice is the dis- 
crete torus A* := ^Z d /27rZ ci . Let w p (x) := e ip ' x where p-x := Y^i=iPi x i- The 
set of functions {w p : p G A*} is an orthogonal basis of 7i (the normalization 
is \\w p \\ 2 = L d ). Therefore the annihilation operators a p = a{w p ) that are 
associated to this basis in the standard way (see [3] or jH]) obey the canonical 
anticommutation relations (CAR) 

{L/ if p — 
otherwise. 

If F is continuous on B := R d /27rZ d , then L~ d E pe A* f (p) ~> f B 7& F (p) 
as L — ► oo. Since we are ultimately interested in the limit L — >• oo, we will 
use the continuum notation even for the finite sums, i.e. we write A, dp F(p) 
for L- d Z peA *F(p),etc. 

Let A be the C* algebra generated by {a+ , a p : p G A*}. For a selfadjoint 
element H <E A, consider the time evolution given by H in the Heisenberg 
picture, i.e. A t := e~ ltH Ae ltH for all A G A. Given a state p on the algebra, 
we define for A G A, Pt(A) := p(A t ). For all t eM, p t is again a state on A, 
and its time evolution is given by the Schrodinger equation 

i£ tPt (A)= Pt ([H,A]). (2) 

We take the Hamiltonian H := H + A$ where 

#o := y dp e(p) a+a p (3) 

1 Thus the fermions are spinless. We could also choose H' A = £ 2 (A 7 C 2 ) as the one- 
particle Hilbert space, to allow for a fermion spin 1/2, without changing the derivation in 
an essential way. 
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is the kinetic energy and 

$ := J dfci . . . dk 4 (h k 2 \®\k 3 k^a^a^a^a^ (4) 

is the interaction. The coefficient function (k\ k 2 | $ | k 3 k 4 ) is antisymmetric 
under exchange of k\ and k 2 and under exchange of k 3 and k 4 , and it contains 
the momentum conservation delta function. For an interaction generated by 
a two-body potential v(x — y) we have 

(h k 2 | $ | k 3 k 4 ) = 6(h + k 2 , h + h) (5) 

- (v(ki - h) - v(k 2 - k 4 ) - v(ki - k 3 ) + v(k 2 - k 3 )) . 



We assume that {k 4 k 3 \ $ | k 2 ki) — (ki k 2 \ $ | k 3 k 4 ). Then $ = $ + . In 
terms of v this condition means that v is real. For simplicity, we shall assume 
that v is symmetric, i.e., 

v(x) = v(—x) . (6) 

We shall call a polynomial F in the creation and annihilation operators 
quartic if it is homogeneous of degree four and contains exactly two creation 
and two annihilation operators. Any such quartic F has a representation 
similar to (0) with a coefficient (ki k 2 | F \ k 3 k 4 ), and we shall always 
assume that the coefficient is given in the properly antisymmetrized form so 
that we can compare coefficients. 

The state pt is determined by its values on monomials in the creation and 
annihilation operators. The two point function in Fourier space is defined by 

Vpq{t) •= Pt(a+a q ). (7) 

We are interested in the Euler scaling limit of the two point function. In 
configuration space it amounts to the rescaling 

x = X/e, t = T/e, e -> 0. 

Recall the Wigner transform of a function ip G L 2 (M. d ) is defined as 

Define the rescaled Wigner distribution as 

X 

e 
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Its Fourier transform in X is given by 



w;&v) = $(v-f)$(y + f). 



We can easily extend this notions to the two point function u pq . In particular, 
we can define the rescaled Wigner distribution Wp(X, V) through its Fourier 
transform: 

Wft,V) = a v+s4 ) . 

Assume that 

W £ (X, V, T) := W s pT/e (X, V) - F(X, V, T) 
as e — > 0. Under the weak coupling scaling assumption, i.e., 

x = X/e, t = T/e, \ = ^/e (8) 
one expects that F(X, V, T) satisfies the nonlinear Boltzmann equation 
dF(X, V, T) 



dt 



+ V-V X F(X,V,T) 



4tt J Ak 2 dhdk A 5 (hi + fc 2 , k 3 + h) 5(E 1 + E 2 - E 3 - E 4 ) (9) 



x \v{k\ — ki) — v{k\ — k 3 



|2 



X 



F kl F k2 (l - F k3 )(l - F kA ) - F k4 F k3 (l - F k2 )(l - F kl ) 



where F kj is short notation for F(X, kj, T) and Ei = e(ki) with k\ = V. If the 
state is homogeneous in space (i.e. translation invariant) at time zero, then 
u pq = 5(p, q)F p (t) for all later times t as well and the Boltzmann equation is 
reduced to 



4vr J dMM^ \v(h - h) - v(h - k 3 )\ 2 (10) 



dF kl 
dT 

xSih + k 2 , k 3 + fc 4 ) 5{E l + E 2 - E 3 - E 4 ) 

x [F kl F k2 (l - F ks )(l - F ki ) - F k4 F ks (l - F k2 ){\ - F kl )\. 



Our goal is to give a heuristic derivation of this equation. 
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Remark 1. Without the symmetry assumption © one has to replace 
the collision kernel \v(k% — k 4 ) — v(ki — k 3 )\ 2 in (jl(Jj) with its symmetrized 
version 

1 2 
-\v(h - k 4 ) - v(k 2 - k 4 ) - v(kx - k 3 ) + v(k 2 - k 3 )\ 

and our derivation remains valid. 

The quartic structure of the collision is due to the quantum nature and 
the weak coupling limit. Instead of the weak coupling limit, one can take the 
low density limit (x = X/e, t = T/e, X = 1 and the density of the particles is 
e). The resulting equation will be the standard nonlinear Boltzmann equa- 
tion where collision term is quadratic with full quantum scattering kernel and 
not just its Born approximation in the weak coupling limit. Technically, the 
emergence of the full quantum scattering kernel can be seen from resumming 
the Born series, see jl] for the simpler case of the Lorentz gas. From our 
experience working on the weak coupling 5J and low density limits in ran- 
dom environments, we believe that a rigorous derivation of the low density 
limit will be somewhat more complicated than in the weak coupling limit. 
However, the key difficulties arising from many-body quantum dynamics are 
already present in the weak coupling limit which we shall focus on. 



3 The equation for the two point function 

The Schrodinger equation can be written as 

(i<9 t - e(p) + e{q))v pq {t) = Xp t (F pg ), (11) 

where F pq := [$, a^a q ] is quartic with 

(h k 2 | F pq | k 3 h) = - 8(q, h){ki k 2 \®\pk 3 ) 

+ 5{q,k 3 )(k 1 k 2 \<5>\pk 4 ) 
+ 5(p,h)(k 2 q\$\k 3 k 4 ) (12) 
- 5(p,k 2 )(k 1 q\<5>\k 3 k 4 ). 

Therefore 

t 

u pq (t) = u^e-^-^ - iX Us e-^- s ^-<^ Ps {F pq ). (13) 



o 
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Thus we need expectation values of quartic monomials, whose evolution 
equation analogous to (|13|) involves ones of degree six, etc. This system 
of equations "hierarchy" ) is similar to the Schwinger-Dyson equations, but 
the commutator structure implies that in an expansion in Feynman graphs 
only connected graphs contribute. 

Let (p ri ,...,r 4 (t) ■= p t (a+a+ 2 a r3 a ri ) and 

Ae(ri, . . . , r 4 ) := e(r{) + e(r 2 ) - e(r 3 ) - e(r 4 ). (14) 

The Schrodinger equation for (p gives 

(xd t - Ae(r x , . . . ,n))(p riy ^(t) = Xp t {[$,a+a+a r ,a r J) (15) 

which integrates to 

<P n ,...,rM = ^^rM^^' 1 -^ 

t 

-iX I d a e-^'^'-"' r ^.([$,a+a+a n , a ». 4 ]). (16) 
o 

Thus 

t 

Pt(F pq ) = Po (F pq ) -ixjds Ps ([$, G pq (t - s)]) (17) 

o 

where G pq = G pq (t — s) is quartic with 
(fa k 2 | G pq | k 3 fa) = e -^-°)^>-M (is) 

(fa k 2 \$\h p)$(fa, q) - (fa k 2 \$\fa p)5(k 3l q) 
-(q k 2 \<S>\k 3 fa)5(fa,p) + (q fa\$\k 3 fa)5{k 2 ,p) . 

Thus the equation (|TT|) for v pq (t) can be written as 

t 

(id t - e(p) + e(q))u pq (t) = Xp (F pq ) - iA 2 J ds p s ([$, G pq (t - s)]). (19) 

o 
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4 Restricted Quasifreeness 



Up to this point, everything was exact; a heuristic derivation of the Boltz- 
mann equation now begins by treating the state p s as quasifree in (fT^j). i.e. 
expressing the term on the right hand side of (fTH|) as a product over u pq (s). 

To this end, it is advantageous to avoid any contractions in the commu- 
tator, i.e. simply leave it in the form [$, G] = $G — G$, which gives 

[<$>,G P g(t-s)) = J dfci . . .dhdh . . ■^Ualal 2 a h a h al,al A a h a h 

M pq (ki, k 2 , h, h, l h l 2 , h, k) (20) 



with 



M pq {k 



1: 



(h k 2 \<S>\h h)(h k 4 \G pq (t - s)\l 2 h) 
-(k 3 h\$\l 2 h)(kx k 2 \G pq {t - s)\U h) 



[21) 



We recall that expectation values of higher order monomials in a quasifree 
state p s can be expressed by the two-point functions (see the Appendix). In 
particular, the four-point function is given by the following determinant 



Vk x h u kil 2 
Vk 2 h V kili 



Ps(atA2 a h a h) ~- 
and the eight-point function appearing in pUj) is 



(22) 



Ps{a ki a k2 a h a h a ks al 4 a l2 a hi 



v kxh v k\l 2 v kxh v kxU 

h> k 2 h v k 2 l 2 h> k 2 h V k 2 U 

U k 3 h V k-ih V k , ih V kzU 

Vk 4 h v ki i 2 v ki i 3 l>k A U 



(23) 



Here each v k i stands for v k i{s) and v ki = —5(k,l) + v k i(s) appears in the 
lower right block because the monomial is not normal ordered. We shall call 
a state p s restricted quasifree if both (J22j) and (J23|) are satisfied. We shall 
assume this condition in the limit A — > 0. 

Return to the derivation of the Boltzmann equation. A Laplace expansion 
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of the determinant gives 



+ 



+ 









Mi i 


l/l T 


fC2tl 


Kill 




VbA In 

K4(3 


Vk A l A 




v k x l 2 




Vk 2 l z 


v k 2 l A 


v k z l x 


V k z l 2 




Vk A h 






"kih 




v k 2 lz 


v k 2 U 




v k A l 2 




"kzh 




v k 2 h 


v k 2 l 2 






Vk x h 


V k z l x 


v k z l 2 




Vk A h 




v k 2 l x 


v k 2 l 2 








Vk 4 h 


v k A l 2 








V k 3 h 


v k 3 l 2 




Vk x h 






v Ul 2 




v k 2 lz 


v k 2 U 



(24) 



+ 



Noting that M pq (ki, . . . , U) is antisymmetric under exchange of l\ and I2 
and under exchange of I3 with U, and that the same is true for each of 
the six summands in the Laplace expansion, we see that we may replace 
every 2x2 determinant by the product of the diagonal elements if we in- 
clude a symmetry factor 4. Moreover, M pg (ki, . . . , / 4 ) is antisymmetric under 
(ki, kit U, I3) — > (&3, &4, I2, h) (this is just the antisymmetry of the commuta- 
tor in its two arguments), but the last of the six summands, ^kzh^k^h^k^^^h, 
is symmetric, so it cancels out. Graphically, this is the cancellation of the 
disconnected term. Thus the term multiplying M pq (ki, . . . , / 4 ) is 



{ v kxlx v k 2 l 2 v k 3 l 3 v k A U 

— Vk^Vkzh^k^Vk^ 

+ ^kih v k A l 2 v k 2 l 3 ^k 3 U 

+ ^h^kzh^kilz^kih 

— v k 2 lx v k A l 2 v k^kzl^) 
{ ^kxlx v k 2 l 2 ^k 3 l 3 ^k A U 

+ {Vk 3 l 2 Vk A l A — Vk 4 l 2 Vk 3 h) { V k 1 h^k 2 h ~ VkxhVfolz)} ■ 



(25) 



Again, the last factor is antisymmetric with respect to an exchange of k\ and 
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k 2 and an exchange of k 3 and k^, so there is another symmetry factor 4, and 
Q 8 ([$,Gpq(t- s)]) = J dk 1 ...dhdl 1 ...dUM pq (k 1 ,...,U) 

^(^xhVkahVkah^U+^kihVkala^kihhau)' (26) 

We remark that if we had used the commutator contraction to express 
G] in ([2 Up. then we would have needed to evaluate only monomials of 
degree six on the state g s . The calculation would have been longer because 
certain cancellations would be less transparent. However this approach has 
the advantage that it requires the quasifree factorization property of g s only 
for degree six monomials instead of degree eight. 



5 Spatial homogeneity 

If we assume that the distribution is homogeneous in space (i.e. translation 
invariant) at time zero, then v vq = 5(p, q)f p {t) for all later times t as well by 
the translation invariance of H. In this case there are further simplifications: 
the term p (F pp ) vanishes and the e(p) — e(q) term in the differential equation 
also drops out. Moreover, for p = q 

(ki k 2 | G pq | k 3 k A ) = e-^-^^-^ihk^lhh) (27) 

x(5(p, h) +S(p, k 3 ) - 5(p, k 2 ) - 5{p, fci)) 

so we can compute the contribution of the first term on the right hand side 
of (ED as 

J dki . . . dhM pq {kx, h)kvkihVk 2 i 2 i>k3hVkdA 

= J dfci . . . dfc 4 e- i(t - s)Ae(fel '-' fc4) (28) 
x8\(h k 2 \$\k 3 k 4 }\ 2 (5(h,p)-5(k hP )) 

X {f hi fk 2 fkz fk A fki fk 3 fk 2 fki ) 

with f p = l- f p . The second term 16M pq (ki, . . . , k^k^k^kAhhah in © 
drops out because with the assignment of momenta it is equal to 

32(«J(A; 1 ,p) - 5(k 3 ,p))(k 1 k 2 \$\k 3 k 2 )(k 3 k i \<5>\hk 1 ) cos[(t - s)(e(k 3 ) - e{ki))], 
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and this quantity vanishes because the delta functions in both $ factors are 
S(ki, k 3 ) and thus S(ki,p) — 8(k 3 , p) = 0. 

Inserting (jSJ) into (|28p. recalling that in our finite volume, 5{p,p) = L d , so 
that S(p, q) 2 = L d S(p,q), and using © and symmetry arguments as above, 
we get 

t 

d t f p (t) = -A 2 J ds Jdk 1 ... dh 8{k t + k 2 , k 3 + fc^e-Kt-'JAeCk!,...,**) 
o 

x2(5(k 4 ,p) - Sih^vik! - h) - v(ki - k 3 )\ 2 (29) 

X(fk 1 ( s )fk 2 ( s )fk 3 (s)fk 4 (s) ~ fk 4 (s)fk 3 (s)fk 2 (s)f kl (s)). 

6 Local approximation in time 

We rewrite the equation as 

oo t 

-\- 2 d t f p (t) = J dE J dse- iE ^ (3(E,p,s) (30) 



-oo 



with 

P(E,p, s) = J dh... dk A 5{k 1 + k 2 , k 3 + h) 2(S(k 4 ,p) - 8{k x ,p)) 

x\v{k x -h) -v(ki - k 3 )\ 2 5(E- Ae(Jfei, . . . , Jfe 4 )) (31) 
x(fki(s)fk 2 (s)fk 3 (s)f k4 (s) - fk 4 (s)fk 3 (s)fk 2 (s)f kl (s)), 

Since v is symmetric (JHJ), f3 is a symmetric function of -E 1 . 

Notice that / and (3 are A dependent and we shall denote them by f x and 
(3 X . We now assume that the limits 

lim/ A (T/A 2 ) = F(T,p), \im 8\E,p,T/X 2 ) = B(E,p,T) (32) 

A — >U A — >{j 

exist and the relation (|3 1|) continues to hold in the limit. We can take the 
limit A — ► in (|3*0*|l and this yields 

oo T 

-d T F(T,p) = lim J dE j ^ e- i£ ( T - s )/ A2 /3 A (£,p, S/A 2 ) . (33) 



-oo 
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We now assume that we can replace the function (5 X by its limit B. Thus we 
have 



T 



-d T F(T,p) = lim j dE J ^ e -iE(T-syx* B ( EiPi S /\ 2 ) . (34) 

-oo 

Interchanging the the integration and performing the E integration, we have 

T 

-d T F(T,p) = hmj ^ B{(T-S)/X 2 ,p,S/X 2 ) . 
o 

Let u = (T — S)/\ 2 . We can rewrite the last equation as 

T 

-d T F(T,p) = lim J du B(u,p,T + u\ 2 ) . 
o 

In the limit A — ^ 0, the right side converges to 

oo oo 

J du B(u,p,T) = l - J du B(u,p,T) = ttB(0,p,T) . 

-oo 

where we have used the symmetry of (3 in E. Combining the last two equa- 
tions, we have derived the Boltzmann equation (jlOl) . 

In this derivation, we used the restricted quasifreeness assumption, spatial 
homogeneity and the existence of the limit for the the two point function v pq 
(cf: (|32j) ). We have not made precise the meaning of the limit and we have 
freely interchanged limits with differentiations and integrations etc. This 
suggests that for a rigorous proof the two point function has to be controlled 
precisely, perhaps through some expansion method. 

The spatial homogeneity of the initial state can be replaced with the 
assumption that the point function at any time scales as 

Our derivation above can easily be extended to this case to give the spatially 
inhomogeneous Boltzmann equation ®. 
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A Quasifree states and determinants 



For finite L, the observable algebra is finite-dimensional, so a state p is 
quasifree if and only if it is given by density matrix coming from a quadratic 
Hamiltonian (see, e.g. 0). That is, 

p(A) = i tr (e-"M) , (35) 

where Z = tr e~ H °. We restrict to states which are invariant under the 
transformations a p — > e ia a p for all a6l For this case we prove below that 
the expectation value of any normal ordered monomial can be computed with 
the following formula: 

/ m m' \ 

P IR II V = ^(-l) m(m ~ 1)/2 det (36) 

\n=l n'=l / - > - 

To simplify notation, we enumerate our finite set A* in some way so that 
we can replace the subscript p G A* by a number i e {1, . . . , A^}, A" = L d . 
Moreover, because (|3T)|) is homogeneous, we may rescale the creation and 
annihilation operators by L _d//2 , so that they obey the CAR a i + a^a i = 
5ij with Sy the Kronecker delta. With these conventions, and by the just 
stated U(l) invariance, 

H = ^2 a tQij a j ■ (37) 

Positivity of p requires H to be hermitian, so = Q^. Thus there is U 6 
U(N) such that Q = UEU* with E = diag {E 1 , . . .,E N }. The operators 
b k = J2j Ukjdj have canonical anticommutation relations b k bf + bfb k = 5m, 
so that 7ik = b k b k satisfies n\ = rik and n^rii = riiUk- Thus e~ H ° is the 
product of commuting factors 

e""° = H (1 + ( e -* - l)n k ) , (38) 
k 

hence Z = H k (l + e~ Ek ), and p(6^fy ) = 5 k i(l + e Ek )~ l . This implies 

piafa,) = U n —^rU a = (1 +e Q Y\, l . (39) 

1 T G 1 
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Because H is diagonal when expressed in terms of b + and b , 



/ m ml 



) 



MIR IK 



(40) 



\fe=i k'=i 



vanishes unless m! = m and (v±, . . . , v rn ) is a permutation of (u\, . . . , u m ) : 
^fc = w 7r (fc). In that case, by the CAR, (|3Uj) equals 



Eq. (J36j) now follows straightforwardly by expressing the product of a + and 
a in terms of b + and 6 and using the definition of the determinant. 

In (J22j) and (|23|). the indices of the annihilation operators are ordered 
downwards, so the factor (— l) m ( m_1 )/ 2 is absent. The procedure of com- 
muting a monomial that is not normal ordered to its normal ordered form 
corresponds to successive row expansions of the determinant in (|23*j) . 
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